Superfluidity of isotopically doped parahydrogen clusters 
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It is shown by computer simulations that superfluid para-hydrogen clusters of more than 22 
molecules can be turned insulating and "solidlike" by the replacement of as few as one or two 
molecules, with ones of the heavier or£/io-deuterium isotope. A much smaller effect is observed with 
substitutional ort/io-hydrogen. Substitutional or£/io-deuterium molecules prevalently sit in the inner 
part of the cluster, whereas ori/io-hydrogen impurities reside primarily in the outer shell, near the 
surface. Implications on the superfluidity of pure para-hydrogen clusters are discussed. 

PACS numbers; 36.40.-c, 67.90.+Z, 61.25.Em 



Clusters of molecular hydrogen are the subject of 
intense experimental investigation, especially after the 
claim of observation of superfluidity (SF) in clusters of 
Af=14-16 P-H2 molecules surrounding a linear carbonyl 
sulfide (OCS) impurity 1]. Experimental data for pris- 
tine or isotopically mixed clusters are not yet available, 
but novel techniques based on Raman spectroscopy may 
soon be able to characterize their superfluid properties 
i- 

Quantitative information and physical insight in the 
properties of clusters of para-hydrogen (p-Ho) molecules, 
both pristine |, i, i, H, 0, i or doped 0, Efl 111 E d , 



has been provided by microscopic calculations based on 
Quantum Monte Carlo (QMC) simulations, making use 
of accurate intermolecular potentials 3,15,3]. In par- 
ticular, a significant superfluid response has been pre- 
dicted at T < 1 K for pristine clusters comprising as 
many as 27 P-H2 molecules [1, [3|- 

To an excellent approximation, P-H2 molecules, as well 
as those of the heavier ort/io-deuterium isotope (0-D2) 
can be regarded as point particles of spin S=0, 0-D2 hav- 
ing twice the mass of P-H2 . A third isotope exists, namely 
ort/io- hydrogen (0-H2), which has the same mass of p- 
H2, but different spin (S=l). Pair-wise interactions be- 
tween molecules of the different isotopes are very nearly 
the same [13]; thus, a mixed cluster of these components 
provides a simple experimental realization of a mixture 
of isotopic bosons. A number of fundamental questions 
can therefore be addressed by studying these systems, 
e.g., how does the presence of one component affect the 
superfluid properties of the other (s), or the structure of 
the cluster ? 

Classical binary Lennard-Jones (LJ) clusters have been 
extensively investigated; for example, it is known that, 
for equal concentrations of the two species, the one with 
the smaller depth of the pair potential well e^j and/or 
with higher particle effective radius cr^j will sit preva- 
lently at the cluster surface [3) [3] ■ Because in a mix- 
ture of, e.g., P-H2 and 0-D2, clj and ctlj are essentially 
the same for all pair-wise interactions, the physics of the 
system is dominated by quantum effects, at low temper- 
ature. 



Path Integral Monte Carlo simulations of such mixed 
clusters have been carried out down to T=2.5 K, but 
without including quantum exchanges [20| . No theoreti- 
cal predictions are presently available for the superfluid 
behavior of this system. Clearly, the sensitivity of the 
superfluid response of a finite cluster to its detailed iso- 
topic composition, has a direct relevance to the design 
and interpretation of ongoing and future experiments. 

In this paper, we present results of a microscopic study 
of P-H2 clusters doped with isotopic impurities; specifi- 
cally, we investigate the effect of the substitution of few 
{Njy < 4) P-H2 with 0-D2 (or, 0-H2) molecules, in {p- 
112) N clusters that are entirely (or, almost entirely) su- 
perfluid. For definiteness, we focus on clusters with N > 
16, at the two temperatures T=0.5 and 1 K. 

The substitution of P-II2 with 0-D2 molecules, gen- 
erally causes pristine superfluid P-H2 clusters to turn 
progressively solidlike, their structure increasingly mim- 
icking that of pure 0-D2 clusters with the same num- 
bers of molecules. For clusters comprising more than 20 
molecules, however, the change from liquid to solidlike 
can occur abruptly. For instance, while the superfluid 
fraction of the pristine (p-H2)25 cluster at T=l K is ap- 
proximately 75%, substitution of a single P-H2 molecule 
with an 0-D2 suppresses SF almost completely. Like- 
wise, two such substitutional 0-D2 impurities suffice to 
suppress the superfluid response (close to 100%) of the 
same cluster at r=0.5 K. 

A much smaller reduction of the superfluid response is 
observed if substitutional impurities are 0-H2 molecules. 
For example, (p-H2)24-(o-H2)i is largely liquidlike at 
T=l K. A crucial observation is that few substitutional o- 
D2 molecules sit in the inner part of the cluster, whereas 
0-H2 impurities are primarily located near its surface. 
These findings underscore the prominent role played by 
long exchanges, especially involving molecules in different 
shells of the cluster, in stabilizing a liquidlike structure of 
large clusters at low T, in turn allowing for SF to occur. 

We model our system of interest as a collection of Nh 
P-H2 and Nd = N — Nh 0-D2 (0-H2) molecules, regarded 
as Bose particles of spin S=0 (S=l), interacting via the 
Silvera-Goldman (SG) pair potential [3]. We study this 
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FIG. 1: (Color online) Superfluid fraction ps^ of the para- 
hydrogen component in clusters of A'^ molecules with one (cir- 
cles), two (triangles) and four (squares) substitutional ortho- 
deuterium molecules, and with one (stars) or£/io-hydrogen 
molecule. Upper panel shows data at T = 0.5 K, lower at 
T = 1 K. Diamonds show results for the undoped clusters 
(entirely superfluid at T=0.5 K). Error bars are comparable 
to the size of the symbols. 
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FIG. 2: (Color online) Radial density profiles for a pure (p- 
112)25 (lower peak at short r) and a (p-H2)24-(o-D2)i cluster, 
at T—1 K. Profiles are computed with respect to the geomet- 
rical center of the cluster. In the case of the doped cluster, 
no distinction is made between molecules of different types. 
Statistical errors are not shown for clarity; they are of the 
order of 5.0x10"* A'^ 



system by QMC simulations, using the continuous-space 
Worm Algorithm 21, 22|; technical details of our simu- 
lations are the same as in Ref. 0- 

Figure [1] shows the superfluid fraction ps^ of the P-H2 
component, as a function of the cluster size TV at T = 1 
K (lower panel) and T = 0.5 K (upper panel). Different 
symbols refer to a different number of substitutional 0-D2 
or 0-H2 molecules (see figure caption). When a single p- 
H2 molecule is replaced by an 0-D2 one, psh is relatively 
little affected for N < 22, with respect to a pristine P-H2 
cluster, but is depressed substantially in larger clusters, 
particularly at higher T. For example, the substitution 
of one 0-D2 in a (p-H2)25 cluster causes ps„ to drop from 
75% to less than 20% at T=l K, while at T = 0.5 K an 
almost complete suppression of the superfluid fraction is 
observed, in the same cluster, when two molecules 
are replaced by 0-D2. It should also be noted that pris- 
tine P-H2 clusters with more than 22 molecules, generally 
display solidlike behavior even when undoped, albeit in 
different degrees [e.g., (p-H2)23 at T=l K] [i,0|. 

Figure [2] compares the radial density profile (at T—1 
K) of the (largely superfluid) pristine (p-H2)25 cluster, 
with that of the (essentially insulating) cluster with the 
same number of molecules, but with one P-H2 replaced 
by an 0-D2 . Proflles are computed with respect to the ge- 
ometrical center of the cluster [i^. The main structural 
change arising from the substitution of a p-B.2 clearly 
occurs in the center of the cluster, where molecules are 
considerably more localized (hence, the higher peak). In 
the vicinity of the surface of the system, density profiles 
are, in fact, almost identical. 

If the dopant is 0-H2, on the other hand, then the 
structure of the cluster and its superfiuid properties are 




FIG. 3: (Color online) Radial density (multiplied by r^) for 
the single 0-D2 and 0-H2 impurities in the doped (p-H2)24-(o- 
D2)i and (p-H2)24-(o- 112)1 clusters, at a temperature of 1 K. 
Profiles are computed with respect to the geometrical center 
of the cluster. 



much less sensitive to the substitution (see Fig. [T]). For 
example, the corresponding profile for the cluster doped 
with one 0-H2 molecule is indistinguishable from that of 
the pristine cluster, on the scale of Fig. O 

An important observation is that a lone 0-D2 molecule 
sits in the central part of the cluster; this is something 
that we observe for all cluster studied, with up to four o- 
D2 substitutional impurities, in agreement with previous 
work The lighter 0-H2 dopant, conversely, is consid- 
erably more delocalized, and indeed is found prevalently 
in the external part of the cluster, as shown in Fig. [31 

The effect of cluster "crystallization," induced in rela- 
tively large clusters by one or two 0-D2 impurities (and 
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FIG. 4: (Color online) Three-dimensional representations of 
the clusters (p-H2)i6-(o-D2)4 (A) and of (p-H2)20 (B). Darker 
color is used for impurity molecules. 

the ensuing depression of P-H2 SF) can also be observed 
in smaller systems, but a greater number of substitutions 
is needed. Figure H] shows the structures of the two clus- 
ters (p-H2)i6-( 0-02)4 (part A of the figure), and (p-H2)2o 
(part B); these pictures were obtained using the proce- 
dure outlined in Ref. [l3|. The pristine cluster has a 
featureless structure, and is entirely superfluid at T < 1 
K. On the other hand, psa is small in the doped cluster, 
whose solidlike structure is evident, with a central axis 
surrounded by rings of molecules. Two of the four 0-D2 
molecules are placed on the axis, the other two on the 
central ring. 

Numerical studies by other authors [13, [3] had al- 
ready yielded evidence of localization of P-H2 molecules 
around a heavy impurity, rendering small, pristine clus- 
ters [e.g., (^-£[2)13] significantly more rigid and sohdhke. 
In all previous works, however, impurities were consid- 
ered such as CO, or HF, not only significantly heavier, 
but, more importantly, featuring a stronger (more at- 
tractive) interaction with the P-H2 molecules, than that 
between the molecules themselves. 

The results presented here, offer insight in the micro- 
scopic mechanism of SF in quantum clusters. In or- 
der for SF to occur, clusters (either doped or pristine) 
must be essentially liquidlike in character, i.e., molecules 
must enjoy a high degree of mobility and delocalization. 
The most important structural difference between clus- 
ters that are insulating or superfluid at T < IK, also 
based on the results for pristine clusters [170|i is that 
the former feature a rigid, solidlike core, possibly with 
some loosely bound molecules on the surface; on the other 
hand, SF is enhanced in clusters whose inner region is 
floppy, with exchanges taking place between molecules 
in the inner and outer shells. 

This is consistent with the notion of quantum melt- 
ing of large clusters, at low T 0, 0], originating from 
permutational exchanges involving all molecules, includ- 
ing those located in the inner part of the cluster. These 
exchanges become increasingly important at low temper- 
ature, where clusters become increasingly liquidlike, and 
consequently superfluid. 

Figure O displays the frequency with which exchange 
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FIG. 5: (Color onhne) Frequency of occurrence of exchange 
cycles of length M (i.e., involving 1 < M < Nh 
molecules) at T=l K. Spikes from left to right refer to a pris- 
tine (p-H2)25 cluster and to the mixtures (p-H2)24-(o-H2)i 
and (p-Il2)24-(o-D2)i . 

cycles of varying length (i.e., involving a different num- 
ber 1 < M < Nh of P-H2) occur (spikes from left to 
right) in a pristine (p-H2)25 cluster as well as in the mix- 
tures (p-H2)24-(o-H2)i and {p-)i2)2i-^ {0-^2)1 at r = 1 K. 
All exchange cycles are clearly suppressed in the cluster 
doped with a single 0-D2 , but the reduction is most dra- 
matic for very long cycles. Conversely, when the cluster 
is doped with an 0-H2, exchanges (other than the very 
longest one) take place at almost the same frequency as 
in the pristine cluster. 

Thus, a single substitutional 0-D2, which is located 
in the center of the cluster owing to its greater mass, 
has a strong inhibiting effect on long exchanges of P-H2 
molecules. As a result, the cluster turns solidlike, and SF 
is altogether depressed. Conversely, a single 0-H2 dopant 
molecule can effectively get out of the way, thereby allow- 
ing for a greater occurrence of long exchanges, including 
those involving P-II2 molecules in the inner and the outer 
shell of the cluster. As a result, the doped (p-H2)24-(o- 
H2)i cluster remains largely liquidlike, with a value of 
ps„ close to 60 %. 

The suggestion that SF is underlain, or strongly en- 
hanced by the occurrence of exchange cycles involving 
also molecules located in the central region of the clus- 
ter is not (as one may naively think) in contradiction 
with the observation that the superfluid density is largest 
at the surface For, as discussed above the primary 
mechanism by which these exchanges promote SF is not 
by locally increasing the value of ps^ , but by stabilizing 
an overall liquidlike phase of the whole cluster. 

Summarizing, we have studied mixed hydrogen clus- 
ters of various sizes, down to T=0.5 K. To our knowl- 
edge, this is the first study of a mixed isotopic bosonic 
cluster including all quantum-mechanical effects, namely 
zero-point motion and permutations of identical parti- 
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cles. Gaining understanding of the effect of substitutional 
impurities, also affords insight into the microscopic origin 
of SF in pristine clusters of molecular hydrogen. 

We find that the supcrfluid fraction ps^ of pure P-H2 
clusters is depressed by replacement of few p-B.2 by 0- 
D2 molecules. The depression is most dramatic for clus- 
ters with more than 22 molcciilcs, which display incipi- 
ent solidlike behavior even when pristine. The reduction 
of the superfluid fraction is considerably smaller in the 
presence of substitutional 0-H2 impurities. Lighter impu- 
rities, which are delocalized throughout the system, have 
less disruptive an effect on long exchanges than the heav- 
ier impurities, which sit in the inner part of the cluster. 
Our findings are consistent with the notion of melting of 
large clusters at low temperature as arising from purely 
quantum-mechanical effects, namely long exchanges of 
indistinguishable particles. 
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